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Abstract 

Local structure of NdFeAsOi-^F^. (x=0.0, 0.05, 0.15 and 0.18) high temperature iron pnictide 
superconductor system is studied using arsenic K-edge extended x-ray absorption fine structure 
measurements as a function of temperature. Fe-As bondlength shows only a weak temperature 
and F-substitution dependence, consistent with the strong covalent nature of this bond. The 
temperature dependence of the mean-square relative-displacements of the Fe-As bondlength are 
well described by the correlated-Einstein model for all the samples, but with different Einstein- 
temperatures for the superconducting and non-superconducting samples. The results indicate 
distinct local Fe-As lattice dynamics in the superconducting and non-superconducting iron-pnictide 
systems. 
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I. INTRODUCTION 



The role of structural topology in the magnetic and superconducting properties of the 
newly discovered iron-based superconductors (FeSC) is one of the active themes of investi- 
gation l|-[3) focussing on finding the mechanisms of the unconventional superconductivity. 
Several studies have shown unprecedented sensitivity of Fermi surface topology to the anion 
(pnictogen or chalcogen) height above the Fe-plane in these materials. The anion height is 
seen to influence the density of states near the Fermi energy and influence the electron-pairing 



properties of the FeSC 



4j-|l0j]. In addition, local structural studies have clearly pointed out 



the importance of the local order-disorder in the spacer layer in the FeSC 
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121. The 



structurally simplest systems among the FeSC, the FeTei_ x Se x ., is found to have lower lo- 
cal structural symmetry than the average crystallographic symmetry, with the Se and Te 
occupying distinct sites Q]- This local inhomogenity is found to have a direct conse- 
quence on the electronic properties of the system |l5j. All these point to the importance of 
the diverging local structure from the average one, putting the FeSC on the same platform 



of cuprates, revealing lattice anomalies and local inhomogeneities [16J, closely related with 
the superconducting properties. Some recent EXAFS studies on F-doped LaFeAsO and 
SmFeAsO have shown a weak anomaly in the Fe-As bondlength fluctuations close to the 



superconducting transition 
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18] . However, such anomalies are found to be much weaker, 



compared to those observed in the cuprates 161 ]. 

Here, we have used extended x-ray absorption fine structure (EXAFS), a fast and 



site-specific experimental tool [19|, |20[], to probe systematically the local structure of 



NdFeAsOi-^Fj. oxypnictides as a function of temperature and charge density. Arsenic 
K-edge EXAFS has been used to retrieve direct information on the Fe-As bondlengths 
in superconducting (x=0.15 and 0.18) and non-superconducting (x=0.0 and 0.05) samples. 
The bondlengths show weak temperature and F-doping dependence, consistent with the 
strong covalent nature of the Fe-As bonds. The corresponding mean square relative displace- 
ments (MSRDs), well described by the correlated Einstein model, reveal different Einstein- 
temperatures (6^) for the superconducting (6s=348±12 K) and non-superconducting 
(B^=326±12 K) samples. The superconducting sample with maximum T c (x = 0.18) ap- 
pears to show a weak temperature dependent anomaly in the Fe-As MSRDs. However, the 
anomaly is well within the experimental uncertainties, similar to case of SmFeAsOo.93iF .o69- 
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These results indicate that there exists, distinct lattice dynamics in the superconducting and 
non-superconducting systems, albeit the differences are smaller. 



II. EXPERIMENTAL METHODS 

Temperature dependent x-ray absorption measurements were performed in transmission 
mode on powder samples of NdFeAsOi_xF x (x=0.0, 0.05, 0.15 and 0.18; the x-values given 
corresponds to the nominal composition, the actual values are lower than the nominal values 



for higher x 21J), at the BM26A beamline 22j of the European synchrotron radiation 
facility, Grenoble (France). The synchrotron light emitted by a bending magnetic source 
was monochromatized by a double crystal Si(lll) monochromator. The samples, synthesized 
by conventional solid state reaction method [2l|], were characterized for their structural, 
magnetic and superconducting properties prior to the experimental run. While the samples 
x=0.0 and 0.05 are non-superconducting, showing tetragonal to orthorhombic structural 
phase transition, the samples with x=0.15 and 0.18 are superconducting with T c =~20 K 



and ~40 K respectively, without any evidence of the structural phase transition [2lJ. For 
the temperature dependent measurements (15-300 K), a continuous flow He cryostat with 
a temperature control within an accuracy of ±1 K was used. A minimum of two scans 
(with high signal to noise ratio) were measured at each temperature to make sure the 
reproducibility. Standard procedure was used to extract the EXAFS oscillations from the 



absorption spectrum 
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20|. 



III. RESULTS AND DISCUSSIONS 

Figure 1 shows Fourier transform (FT) magnitudes of the arsenic A"-edge EXAFS oscilla- 
tions (the insets show the EXAFS oscillations), measured on NdFeAsOi_ x F x . (for x=0.0, 
0.05, 0.15 and 0.18) samples, providing partial atomic distribution around the arsenic 
atoms. At room temperature, all the samples studied have tetragonal structure (space 



group P4/nmm) 21j. The x=0.0 and 0.05 samples show structural phase transition from 



tetragonal to orthorhombic phase below ~ 130 K [2lJ. There are four Fe atoms as the 
near neighbors of arsenic at a distance ~ 2.39 A (main peak at ~2 A). The next nearest 
neighbors are four Nd atoms at a distance ~ 3.28 A and four oxygen atoms at a distance ~ 
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FIG. 1: Fourier transform (FT) magnitudes of the arsenic K-edge EXAFS oscillations for the 
NdFeAsOi-^F-c (x=0.0, 0.05, 0.15 and 0.18). Corresponding EXAFS oscillations are shown in the 
insets. The FTs are not corrected for the phase shifts, thus represent raw experimental data. A 
systematic temperature dependence is evident from the gradual change in the FT intensities. 
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FIG. 2: Fourier transforms of the arsenic K-edge EXAFS oscillations measured (symbols) on 
NdFeAsOi_ x F x (x=0.0, 0.05, 0.15 and 0.18) at low temperature (12 K), together with a single 
shell fit (solid line). The filtered EXAFS oscillations corresponding to the first shell (As-Fe pairs) 
together with the model fits (solid lines) for x=0.0 and x=0.18 samples at 12 K are shown in the 
inset. 

3.51 A, followed by eight As atoms at ~ 3.9 A. Contributions of these distant shells appear 
mixed, giving a multiple structured peak at ~2.5-4.5A (Fig. 1). Apparently, the atomic 
distribution around the arsenic is similar in all the samples, as evident from the FT as well 
the EXAFS oscillations. 

In the NdFeAsOi_ x Fa; system, the arsenic atoms have the nearest neighbors as Fe atoms, 
and their contributions to the arsenic K-edge EXAFS is well separated from all other shell 
contributions. This makes the arsenic K-edge EXAFS data highly suitable for extracting 
quantitative information on the Fe-As bondlengths and the related mean square relative 
displacements (MSRDs). Here we have exploited this possibility and performed single shell 
model fits to the EXAFS oscillations due to the Fe-As bondlengths. 

The EXAFS amplitude depends on several factors and is given by the following general 
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equation [19|, |20( : 

i 1 

where Nj is the number of neighboring atoms at a distance Rj, S 2 , is the passive electrons 
reduction factor, fj(k,Rj) is the backscattering amplitude, A is the photoelectron mean free 
path, 5i is the phase shift, and of is the correlated Debye- Waller (DW) factor, measuring 
the mean square relative displacements (MSRDs) of the photoabsorber-backscatterer pairs. 
Apart from these variables, the photoelectron energy origin, Eo, is another input needed 
for the modeling of the EXAFS. In the present case, we have used a single-shell EXAFS- 



modeling 



191 ] to extract the Fe-As bond correlations. For such a single-shell analysis, we have 



used the WINXAS package [23| , with backscattering amplitudes and phase shifts calculated 



using FEFF [24| with the crystal structure data from diffraction 21] as input. In the single- 
shell least-square fit, the number of independent data points, N in d ~(2AkAR)/7r [19| were 
about 13 (Ak=14 A -1 (k=3-17A~ 1 ) and AR=1.5 A), but, we used only two parameters, 
the radial distance R« and the corresponding MSRD of, as the refinable variables, fixing all 
other parameters (Sq=1), to obtain reliable information on the Fe-As bond correlations. The 
uncertainities in the derived two paramters, R and o 2 , were estimated by standard EXAFS 
method [20]. The canonical approch to error estimation is to determine the region around 
the best fit that contains the true value with a certain probability 0. The projection of that 
volume onto an axis corresponding to a parameter gives the paramter errors, for a chosen 
value of (3. In general this uncertainities dependent on the experimental as well as EXAFS 
data extraction procedures, in addition to the uncertainties coming from the statistical x 2 - 
procedure. In the present case, due to the adoption of identical experimental conditions and 
EXAFS data extraction procedures, the first kind of uncertainites are minimum. However, 
we have set the error bars on derived parameters, R and o 2 , to two times the highest 
uncertainity estimated. This is to underline the fact that we are interested in discussing 
the relative changes rather than the absolute values of these paramters. Figure 2 shows the 
single shell model fit in the real and k-space at 12 K for samples with different x. Here 
it should be mentioned that the samples used in the present study are phase pure with 
the amount of impurities below the sensitivity of the x-ray diffraction. The main impurity 
)hase present in the doped samples is NdOF phase, which showed a F-doping dependence 



2l|. However, in the present case, this impurity phase no influence on the arsenic i^-edge 
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FIG. 3: Temperature dependence of the Fe-As bond distance for the NdFeAsOi-^F^ system. The 
Fe-As bond-length show only weak doping and temperature dependence, indicating covalent nature 
of the bond. However, there appears a small anomaly (arrow marks in the upper panels) in the 
superconducting samples close to the transition temperature, with no such changes evidenced in 
the non-superconducting compounds. 



EXAFS data. The impurities which possibly can inter 



ere with the present EXAFS results, 



2l| 



thus ruling out any contribution 



like, the FeAs phases, NdAs phases were less than 1% 
from such impurity phases in the data presented here. 

Figure 3 shows temperature dependence of the Fe-As bond distances determined from 
the single shell analysis of the arsenic K-edge EXAFS data. The Fe-As distance in all the 
samples are very similar, around 2.39 A. The over-all temperature dependence of the distance 
also looks very similar. However, one can notice a small change in the superconducting 
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samples somewhat close to the transition temperature, with no such changes evidenced in 
the parent compound. An earlier temperature dependent x-ray diffraction study on the 
NdFcAsOo.85 superconductor has showed an abrupt change in the Fc-As distance around 
T c {27]]. Present results also show that the Fe-As distance tend to change near the T c for 
the superconducting (x=0.15, 0.18) samples (Fig. 3 upper panels). But, we do not see 
significant changes in the Fe-As bondlengths or corresponding MSRD across the structural 
phase transition for the parent compound, which is consistent with earlier EXAFS studies 



on LaFeAsO 
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281 ]. SmFeAsO [18[ and BaFe2As2 [29(|, showing no change in the Fe-As 
bondlength and corresponding MSRD across the structural phase transition. As shown in 
a recent diffraction study of the REFeAsO (RE stands for rare-earth) system, data from 
single crystals are important in for a better understanding of structural phase transition 
properties 



251 ] . For example, a recent study on the NdFeAsO single crystal revealed two low- 
temperature phase transitions in addition to the tetragonal-to-orthorhombic transition at Tg 
~142 K [26j indicating the importance of single crystals in revealing the intricate properties 
of FeSC. Naturally EXAFS studies using single crystals are also going to be important, in 
addition to the angular dependent local information provided by the polarization dependent 



EXAFS measurements 



16| 



Figure 4 shows mean square relative displacements (a 2 ) of the Fe-As pair, describing 
distance-distance correlation function (correlated Debye Waller factors). The MSRD is a 



sum of temperature independent (<Tq) and temperature dependent terms [l9j, i.e., 

a Fe-As = °0 + a Fe-As(T) 

In several cases, the temperature dependent term can be described by the correlated 
Einstein-model Q,Q|, 

a Fe _ As (T) = - —coth{- 



Z^Fe-As^E 2ksT 

where Fe-As is the reduced mass of the Fe-As bond and u F is the Einstein- frequency. The 
related Einstein-temperature can be obtained from the expression = hidE/hs- In the 
present case, temperature dependence of a Fe _ As is found to follow the correlated-Einstein 
model. Description of a 2 Fe _ As using this model yield Q F values to be 348±12 K and 326±12 
K respectively for the x=0.18 and the £=0.0 samples. Notice that the uncertainity in 
given above are obtained from the least-square fit considering the error bars on the data 
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FIG. 4: Temperature dependence of the Fe-As MSRD (symbols) for the NdFeAsOi-^F^ (x=0.0, 
0.05, 0.15 and 0.18) system. Upper panel shows the data on the superconducting samples (open and 
filled symbols corresponds to x=0.18 and 0.15 samples respectively), while the lower panel shows 
the data on the non-superconducting samples (open and filled symbols corresponds to x=0.0 and 
0.05 samples respectively). The dotted lines represent the correlated Einstein model fit to the data. 
The Einstein-temperatures, O^:, are 348±12 K and 326=L12 K respectively for the superconducting 
(x=0.18) and non-superconducting (x=0.0) samples. The temperature independent Oq is found to 
have similar value (~0. 00073) in both the cases. 

poitns. The obtained Einstein-frequencies are in the range of phonon- frequencies observed 



in Raman studies for the modes involving arsenic and iron atoms [3l|, |32|. In addition, 
different Einstein frequencies (ooe) indicate different local force constants (k=fiFe-As^E, 
where k is the effective force constant) for the Fe-As bonds in the superconducting and 
non-superconducting samples. The calculated local force constants for the Fe-As bonds are 
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~6.65 eV/A 2 and ~5.85 eV/A 2 respectively for the superconducting and non supercon- 
ducting samples. Thus the Fe-As bondlength seems to get harder in the superconducting 
regime. This observation is an indication of non-negligible role of the lattice modes in the 
superconductivity of these materials. 
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Comparing the already available local structural data on different oxypnictides 17], 



18 



28 



33 



34| one can see that the Fe-As bond in these systems shows onlylittle changes with 



doping and temperature. An EXAFS study on a series of oxypnictides [33] have shown that 
the Fe-As bondlength and the related MSRDs hardly show any change with the varying rare 
earth size, consistent with the strong covalent nature of the Fe-As bonds. On the other hand, 
earlier studies on F doped La-1111 [n| and Sm-1111 [3] compounds have shown presence 
of an anomaly in the temperature dependence of the Fe-As MSRDs. But these anomalies 
were seen to be very weak comparing to that observed for the cuprates jl6]. Indeed another 



study on superconducting La-1111 system 



did not indicate any such anomaly. Such 



29] and F doped Ce-1111 



anomalies were also not seen in the K doped Ba-122 system 
system 30J. From the results presented in Figure 4, for the F-doped Nd-1111 system, there 
is no evidence for a significant anomaly in the MSRD of Fe-As bonds associated with the 
superconducting transition. Here, we make an explicit comparison of the Fe-As MSRDs 
of the superconducting Sm-1111 and Nd-1111 samples to underline that the Fe-As MSRDs 
tend to change around the superconducting transition temperature, albeit the anomalies are 
indeed weak in compare to what has been seen in the copper oxide superconductors [tgI ]. 

Figure 5 shows temperature dependence of the MSRDs of the Fe-As bond for the Nd- 
FeAsO, compared with that of the LaFeAsO 28 1, SmFeAsO Q and BaFe 2 As 2 29]. In 
all the cases, the MSRD data are extracted from the arsenic .fT-edge EXAFS. Correlated- 
Einstein model fit to the data is also indicated in the figure (as dotted lines). From Fig. 
5, it can be seen that the Fe-As MSRDs show almost identical temperature dependence for 
these different oxypnictide systems. Although the temperature dependence looks very sim- 
ilar, there seems to be a systematic change in the Einstein-temperature. The 0# values are 
respectively, 316±5 K, 318±10 K, 326± 12 K, and 328± 12 K for the LaFeAsO, BaFe 2 As 2 , 
NdFeAsO and SmFeAsO. There is a clear increase in the force-constant of the Fe-As bonds, 
with increasing RE size in the "1111" series. Optical studies using single crystals of the 
"1111" series reveal that the Fe-As stretching mode in these oxypnictides shows substantial 
hardening with the rare-earth size, being harder for the NdFeAsO and SmFeAsO showing 
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FIG. 5: Temperature dependence of the Fe-As MSRDs for the NdFeAsO system (open circles in 



all panels) is compared with that of the LaFeAsO 



28( (filled diamonds (a)), SmFeAsO [l^ (filled 



triangles, (b)) and the BaFe2As2 [29] (filled circles, (c)). For facilitating a better comparison, the 



LaFeAsO 



2$, BaFe 2 As 2 |29| and SmFeAsO [18[ data are scaled respectively by 0.0004, 0.0008 



and 0.0017. Panel (d) presents the temperature dependence of the Fe-As MSRD for the F doped 



superconducting Nd-1111 and Sm-1111 



181 ] (with the latter scaled by 0.0012). Overall temperature 



dependence of the two superconducting systems are identical. Inset in panel (d) shows a zoom 
around the superconducting transition temperature. The error-bars in the inset are identical as 
that given in Ref. 
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higher T c in compare to the LaFeAsO [35]. The E u mode frequency for the BaFe2As2 



is 



also found to be similar to that of the "1111" series 



35]. In the "1111" family, the changing 



rare-earth ionic size leads to a change in the pni ctog en height above the Fe-plane 



thus a different interlayer atomic correlations 



11 



33] and 



12 1 . Local structural studies using atomic 



pair distribution function analysis on Bai_ x K E Fe2As2 show that, although the changes in the 
e-As bonds are minimum, the FeAs4 tetrahedra show a systematic change with K doping 



34J. Indeed a small change in the pnictogen height is known to have significant effect on 



the electronic properties of the system through changing degeneracy between different Fe 
3d bands (in particular between the d x 2_ y 2 and d xz /d yz ), with a direct implication on the 
magnetic structure and superconductivity Q . Very recent photo-emission studies have 
showed this pnictogen height dependent changes in the electronic structure 37] confirming 
the importance of even very small local structural changes in determining the properties of 
the FeSCs. 

To further enlighten possible anomaly around the superconducting transition temper- 
ature, we have made an explicit comparison between the temperature dependent Fe-As 
MSRDs for the F doped highest T c Nd-1111 and Sm-1111 18] superconductors (see, e.g., 
Fig. 5[d]). The inset in figure 5[ d] i s a zoomed view around T c . The error-bars given 



here are identical as given in Ref. 



181 ] . For the NdFeAsO0.s2F0.1s sample, there is a weak 



anomaly in the temperature dependence of the Fe-As MSRDs similar to that observed in 
the superconducting SmFeAsO .93iF .069 yM- However, the anomaly is well within the ex- 



perimental uncertainties. This situation is different from that of cuprates 16] where several 
measurements showed clear lattice anomalies associated with the superconducting transi- 
tion. Although more experiments are needed, the weakness of the anomalies does not imply 
that lattice fluctuations have minor role in the Fe-based superconductors. Indeed, the an 
empirical relationship between the superconducting transition temperature and the FeAs4 



tetrahedra in FeSCs 



re superconducting 



341 ] . The increased 



4j-[6j clearly establishes the role of local lattice in t 
properties. Local structural studies also underline the same point [33, 
force constant of the Fe-As bonds in the superconducting sample compared to the parent 
compound in the Nd-1111 system indicate the importance of lattice dynamics in determin- 
ing the properties. Systematic changes in the Fe-As force constant within the different 
iron-pnictide compounds, as revealed by the present study, further support this. 
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IV. CONCLUSIONS 



In conclusion, we have carried out systematic temperature-dependent local structural 
studies of the NdFeAsOi_ x F z (x—0.0, 0.05, 0.15 and 0.18) high temperature pnictide super- 
conductor using arsenic K-edge EXAFS. In all the samples, the temperature dependence of 
the mean-square relative-displacements (MSRDs) of the Fe-As bonds are found to follow the 
correlated-Einstein model. However, the Einstein-frequency for the superconducting sample 
is higher than that of the parent compound, indicating a hardening of the Fe-As bond in the 
former compared to the latter. The over-all temperature dependence of the MSRDs of the 
Fe-As bond seems to be similar in the NdFeAsO, SmFeAsO, LaFeAsO and BaFe2As2, but 
with a systematic variation of the corresponding force constants. The changes occurring in 
the FeAs4 tetraheadra, together with the coupling between the active layer and the spacer 
layer, may account for the changes in properties of different iron oxypnictide superconduc- 
tors. 
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